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Abstract 

The rare decays B~ t i>, B ^ t~ ^ X vu and h ^ X all con- 

tain third generation leptons in the final state, and hence are sensitive to new 
physics that couples more strongly to the third family. We present model 
independent expressions for these decays that can be useful to study several 
types of new physics effects. We concentrate on super symmetric models with- 
out R-parity and without lepton number. We also assume a horizontal U{1) 
symmetry with fermion horizontal charges chosen to explain the magnitude 
of fermion masses and quark mixing angles. This allows us to estimate the 
order of magnitude of the new effects, and to derive numerical predictions 
for the various decay rates and for the forward-backward asymmetry and the 
r polarization components measurable m. h ^ X t~~ . In some cases the 
branching ratios are enhanced by more than one order of magnitude, render- 
ing foreseeable their detection at upcoming B-factories. We also discuss how a 
measurement of asymmetries in 6 — > X t+ can be crucial in distinguishing 
between different sources of new physics. 
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I. INTRODUCTION 



The standard model (SM) of the strong and electroweak interactions provides a success- 
ful description of all the phenomena involving the known elementary particles. However, 
experimental results involving fermions of the third generation are far less precise than for 
the first two generations. From the theoretical point of view, several models of new physics 
predict larger deviations from the SM for processes involving third generation fermions [|l|. 
This is also the case in a class of supersymmetric (SUSY) models without R-parity where 
violation of lepton (L) and baryon [B) number couples more strongly to the heavier fermions 

In recent years, the experiments at the CERN e'^e collider LEP have provided us 
with most of the results on b decays into the third generation leptons r and z/,- . This is 
because the LEP environment has the advantage over symmetric 5-factories (like CLEO) 
or hadron colliders (like CDF) of allowing accurate measurements of the missing energy 
associated with primary h ^ v^- oi secondary h ^ t ^ final state neutrinos. In this way, 
decay modes yielding a missing energy spectrum harder than the usual semileptonic decay 
can be effectively measured or constrained. At LEP, measurements of an excess of events 
over the semileptonic background with missing energy between 10 and 30 GeV f^-Q was 
interpreted as the signature of the decay B ^ X^t followed hy t ^ v X . This yielded 
BR(6 -> X^rz/) = 2.68 ± 0.34% in agreement with the SM prediction BR(6 X^rv) = 
2.30 ± 0.25% |jlO|- Using a large missing energy tag the L3 Collaboration set the 90% 



confidence level upper limit on the exclusive leptonic decay B tv |T1|] 

BR(S- ^ r P) < 5.7 X 10"^ . (1.1) 



In ||T2| it was discussed how similar analyses can yield a limit on the flavor changing decay 



h ^ X u u . Based of the full LEP-I data sample, the ALEPH Collaboration derived a 



preliminary 90% confidence level limit on this decay mode [13 



BR(6-^Xz/i?) < 7.7x 10"^ (1.2) 



Being only one order of magnitude above the SM predictions, the limits ( |1.1| ) and ( |1.2|) 
imply strong constraints on several models of new physics |T^Jl^ . 

The tight limits on very large missing energy events (-Emiss > 35 GeV ) in h decays 
reported by the ALEPH collaboration allowed to estimate order of magnitude bounds 



on BR(i?(f ^ r+ r ) and BR(i?5 ^ r+ r ) |T3|. In Ref. [|r^ it was also argued that a weak 
upper limit on the branching ratio for 6 ^ X t+ of the order of the semitauonic branching 



ratio is implied by the LEP missing energy measurements. The limits were estimated as |T3 

BR(5d ^ r+r-) < 1.5 x 10"^ 
BR(fi, ^ T+T-) < 5.0 X 10^2 

BR(6 Xt+ T-) < 5.0 X 10-2 . (1.3) 

Even if several orders of magnitude above the SM rates, these figures still yield bounds on 
some new physics parameters which are unconstrained by other processes [^. 

In the near future, experiments at i?-factories will reach a much higher sensitivity in 
the study of 6 ^ r , i/^ decays. A measurement of B^ ^ rv appears to be accessible even 
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at the low SM rate. With refined experimental techniques and after few years of run, the 
experimentally very challenging decay h ^ X v v could also be measured, at least in some 
exclusive decay channel. Because of the even lower rates and of the expected low efficiency 
in r identification, the decays t+ and b X might be out of the reach of 

experiments like BaBar and BELLE if their rates are at the SM level. However, as we will 
discuss, some new physics models predict decay rates more than one order of magnitude 
above the SM. Therefore, dedicated studies of 6 — > r , i/,- decays at 5-factories represent a 
powerful tool for detecting signals of new physics. 

The paper is organized as follows. In section II we introduce general four fermion ampli- 
tudes for the decays lu , Bg ^ l~ and b ^ XgU u {q = d, s) and we give the results 
for the various branching ratios. Since our approach is essentially model independent, it is 
well suited to study different types of new physics contributions. 

In section III we concentrate on the decay b ^ X l~ . To take into account the effects 
of new physics, the standard basis of operators contributing to the effective Hamiltonian of 
the decay [|16| has to be enlarged. We generahze it by introducing a set of operators for 
the right-handed flavor changing current sn'-^^bn together with a set of new scalar operators. 
These new effective operators arise in several new physics models, like SUSY models without 
R-parity 0, models with leptoquarks |T^, left-right symmetric models |T^, etc. We study 



the various observables measurable in the decay: the inclusive rate, the forward-backward 
asymmetry and the r polarization asymmetries. Since for the decay channel b X 
the effect of the r mass is non negligible and the average energy of the final hadronic 
system is not very large, in our computation we retain all the fermion masses. We next 
apply our results to the study of SUSY models without R-parity and without L number. 
The theoretical framework is presented in section IV. Since in these models the values 
of the various R-parity violating couplings is not determined, without further theoretical 
input no numerical prediction of the corresponding effects on B decays is possible. In 
order to estimate the new physics effects, we appeal to models where the magnitude of the 
fermion masses and CKM mixing angles is explained by assuming some horizontal U{1) 
symmetry. This framework provides us with additional theoretical constraints yielding a 
set of numerical predictions for the various L violating couplings, and allowing for order of 
magnitude estimates of the various decay rates. 

In section V we present a numerical analysis of two representative models and we discuss 
the results. In particular, we compare our estimates for decays involving the transition 6 — >■ r 
with the corresponding decays involving muons, and we confront the predictions of our new 
physics models with the SM. Finally, section VI contains the summary and our conclusions. 
The set of input parameters used in the numerical analysis is collected in an Appendix. 



II. THE DECAYS ^ tu, B ^ t+ AND b^Xvu. 



The most general non-derivative effective four-fermion interaction involving a b quark, a 
q = d ,s or u quark, and a pair of leptons i and i' can be written in the form 

'His = E(^r,6) {i[CaTa +a'r,75]f) (2.1) 

a 



3 



where Ta = {/, 75, 7'^, 7^75, o"^'^} with a = {S, P,V, A,T} the standard basis of operators 
of the Chfford algebra. In (|2.1|) we have factored out the Fermi constant Gp so that all 
the coefficients Ca and C'^ are dimensionless. Even in the presence of new physics, most 
of the rare B decays depend only on a subset of the operators in (^TTI) . This is due to 
the fact that for purely leptonic B decays, several matrix elements of the quark operators 
vanish. Assuming that neutrinos are described by two component left-handed spinor fields, 
the number of relevant operators is further reduced when neutrinos appear in the final state.|^ 
We will now list the general expressions for the different decays. 



A. The decay B 



lu. 



The decay B^ lu is described by the effective Hamiltonian ( PTT| ) with q = u , £ = k 
and i' = Uj where i,j = 1, 2, 3 correspond to the different lepton flavors. In the presence of 
new physics (for example in SUSY models without R-parity) i 7^ j is an open possibility. 
Since flnal state neutrinos are not detected, in these cases a sum has to be taken over all 
the allowed decay modes. 

The general amplitude for this decay involves a set of matrix elements (0| gF^ 6 . They 
vanish for the parity-even operators F5 = J and Fy = 7^ due to the pseudoscalar nature 
of the B meson. The tensor operator Ft = cx^i/ is antisymmetric in the Lorentz indices, 
and hence its matrix element must vanish as well, since the only available four-vector is the 
momentum of the B meson. Therefore, only the matrix elements of the pseudoscalar 
and axial- vector operators contribute. They are given by the PC AC (partial conserved axial 
current) relations 



{0\uYl5b\B-)=tfBp'l 
{0\u-f5b\B-) = -ifB 



nil 



nib + niu 



TJlb 



(2.2) 



Under the assumption of two-component left-handed neutrinos = P^v , with P/, = |(1 — 
75) , we further have I^^^ul = —Jul and J'-y/^'y^UL = —I'J^i'l ■ For on-shell flnal state leptons, 
the latter operator contracted with the B meson four-momentum = + kp yields 
Pb (J'lu^L ) = nil (J^L ) • Hence the amplitude for the B^ I v decay reads 



J^^ = ifj^mB Gf 



rriB 



(2.3) 



The corresponding expression for the decay rate is 



BR{B- ^ W) = flTB 



GlrriBTn 



IGtt 



1 - 



m'i 



E 



{Ca - CT') - —{Cp - Cp') 

Tfli 



(2.4) 



*We assume that even in the presence of new physics, the neutrinos produced in B decays are 
mainly the SM ones. This is not a strong assumption, since light right-handed neutrinos are 
theoretically disfavored. If the 'SM neutrinos' have non-vanishing masses, operators that vanish in 
the massless hmit are suppressed at least as rriy/mb and hence always negligible. 
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where the sum over the index i accounts for possible z/j 7^ ui decay channels. For / = r 
new physics can induce sizeable enhancements over the SM rate (this can occur for example 
in multi Higgs doublet models where new contributions arise from charged Higgs exchange 
diagrams [0). However, a large theoretical uncertainty in predicting the branching ratios is 
associated with the present poor knowledge of fs ■ Therefore, it could be difficult to identify 
unambiguously new physics effects in this decay. 
In the SM = Cj?' = and 



SM 



-72 K 



ub • 



(2.5) 



Using the set of reference parameters listed in the Appendix, for the SM branching ratio we 
find BR^^(5- ^ r z/) = 7.1 X 10"^ . 



B. The decay Bq 



l+l- 



A detailed analysis of the decay Bq and of the possible types of new physics 

contribution was presented in For on shell r's, (^7^0 = (^r+ + ^t-) (^7^*0 = so that 
also the contribution of the axial- vector operator (0| (77^75 6 |-B) vanishes when contracted 
with the leptonic vector current. The general form of the amplitude reads 



= ifsg rriB Gf 



ri ^q 
^ ^A 



rriB 



{ll,l)+Cl' {II) 



(2.6) 



and the corresponding branching ratio is 

BR{Bq 



^ M = JbJb 



;7r 



\ 



1 - 



4m? 



m. 



"2 mi 
rriB 



'^A 







to 













• (2.7) 



In the SM, C|> and Cj, arise from penguin diagrams with physical and unphysical neutral 
scalar exchange, and are suppressed as ~ (m;,/mvi/)^ ||2^. The decay rate is then determined 
by 



CI 



(yV,lVtq 



SM 



^8 TT sin^ 6w 

where Xt = mf/m^^, and at leading order ^1 



(2.8) 



YJx) 



X 



X - A 



+ 



3x 



In X 



X 



[X 



(2.9) 
= 9.1 X 

We stress that as for ^ ru, also for this decay theoretical predictions 
are plagued by the large uncertainty in the value of /^^ , which can easily mask new physics 
effects. 



Using the parameters listed in the Appendix, we find BR {B, 

10-' \Vtq/Vt, " 



r ' T 
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C. The decay b 



Xq V V. 



The decay h ^ XqV v was thoroughly studied in |T2[ . In the presence of new physics, 
the flavor of the two final state neutrinos can differ. Still, under the only assumption that 
neutrinos are purely left-handed and effectively massless, the general form of the amplitude 
has the remarkably simple form [O] 



A 



(2.10) 



In terms of the coefficients in ( |2.1| ) we have Cl,r = [{Cy — C'y)±{C a — C'^)]! A. Summing over 
the undetected neutrino flavors the branching ratio normalized to the semileptonic decay 
reads 



BR(6 



XqVv) 



J2ij ( 


\cf 


2 + 






8 






c/^b 


) 



BR(6 -> X^eu) 



(2.11) 



where fps{x) = 1 — 8a; + 8x'^ — — 12a;^lnx ~ 0.5 for x = ml/ml is the phase space 
factor for the semileptonic decay. In the SM the decay proceeds via W box and Z penguin 
diagrams, and only one operator contributes to the decay: = {qLl^bi) 

The corresponding coefficient reads 



ci 



SM TT sin^ Ow 



where 21 



X 



2 + X 3x — 6 
7 + 7 ttt; In X 



X 



\x 



The additional 1/m^ and as 



corrections to this result can be found in [ll2|,| 



(2.12) 



(2.13) 



In con- 



trast to the previous decays, theoretical predictions for b ^ XgU u are remarkably free from 
uncertainties. In fact all the parameters entering in (|2.11|) and ( p.l2|) are known with good 
accuracy (the main uncertainty comes from rrit), there are no long distance effects and QCD 

From the theoretical point of view, new physics affecting this 



corrections are small 12 



decay could be identified in a very clean way. At leading order, the SM prediction for the 
>SM,, , ^ ^.^-A -4.4 X 10-5 \VtJVtsf. 



branching ratio is BR^^(6 



XgUU) 



III. THE DECAY h^Xl+l-. 

A. General operator basis. 

In the SM, the effective Hamiltonian for the weak decay 6 — /~ is defined in terms 
of a set of ten effective operators Oi-Oio ITHl.f] At leading order, and neglecting small 



^For simplicity we will restrict ourselves to the case when the final hadronic system carries 
strangeness {X = Xg). Generalization to the case X = X^ requires introducing the additional 
operators O" = {dfl'^^bf) (u^JfiU^) and O2 = {dLlfj-^L) (^l^/^'^l) ' keeping the terms proportional 
to Vub and including new long distance effects. 
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contributions induced only through operator mixing, the operator basis can be truncated to 
the following set |24 



Oio = {sl-f.bD {h.l^l) ■ (3.1) 

In the models we want to study, a larger set of non-renormalizable operators arises. As it 
will become clear in the next section, after Fierz transformation squark exchange induces 
at the tree level the new operators O2 , Og , O^q which are analogous to the corresponding 
operators in (|3.1| ) with the replacement (sl7^6l) ~^ {sR'y^bR). While O9 and 0[q contribute 
directly to the decay, O2 enters only at the one loop level, so that the corresponding short 
distance contribution is small. However, O2 induces also new long distance effects associated 
with cc resonances which can further enhance the decay rate above the SM. For this reason 
we include O2 in our set. At the new physics scale m > 100 GeV where these operators are 
generated, other operators from new physics appear only at the loop level. Since they give 
only suppressed short distance contributions we set to zero the corresponding high energy 
coefficients. However, in the evolution from the scale m down to operator mixing occurs, 
and from the point of view of the low energy theory a clear distinction between tree-level 
and loop-level contributions is lost. This forces us to extend the basis to the following set 

o[ = {sii.bD (4^,4) 

O'j = (sV^,[mbPi + m,PR]6")P^^ 

0[o = {sli.bD {h.l^l) . (3.2) 

In addition to the new set {0[} , the exchange of sleptons induces at the tree level new scalar 
operators which also contribute directly to the decay 

{slbD ill) 

(slbt) (Hi) 
{sWr) ill) 

{sWr) ih,l) . (3.3) 

These operators do not mix with {Oi} and {O^} , neither with new loop-induced four-quarks 
scalar operators, which at lowest order do not contribute to the decay. Hence we can truncate 
the basis to the subset of operators listed in ( p.l|) , (|3.2|) and ( p.3|) . 

Including the new physics, we can now write the Hamiltonian density as 

Kfr = -Gp\{s^^{C,PL + C',Pn)b) (l^l) + {s^,{C,oPl + C[,Pn)b) (l^^H) 



o 
o 



s' 

9 

s ' 

10 



7 



2i'^[sa^,[ {CrPn + C^Pi)m, + (C^Pl + C'.Pr) 



+ ( rs {CIPl + C'iPR)h) ill) + ( s {C!,Pl + C'ZPR)h) {Hi 



(3.4) 



where we have factored out the Fermi constant Gp so that the coefficients C are dimension- 
less. 



The SM contributions to (|3.4|) enter through the coefficients 



C, 



7,10) 



and 



TT 



V V* p"^^ 
^tb ^ts , 



(3.5) 



where Cj are the usual QCD improved Wilson coefficients [^. In the leading logarithmic 
approximation, and neglecting small operator mixings, we have P^-p^] 



23 — 23 



_6_ _ 12 

fj23 -\- fj 23 



Ci[mb 

C2{mb 



C2{Mw) 
C2{Mw) 



58 , _io 29 , 

(ri 23 - 1) [ri- 

135^' ' 189^' 



28 
23 



An 



4 , 11, 8 , 



C2(Mh.) 
C2{Mw) 



cioirrib) = Cio(M, 



(3.6) 



where rj = as{Mw)/c(s{mb) and at the renormalization point /i ~ Mw , ci{Mw) = and 
C2{Mw) = 1- The remaining three coefficients are functions of rrit/Mw and their explicit 
expressions can be found in |2^. Using rrit = 176 GeV and sin^ 9 = 0.23 (where 9 is the weak 
mixing angle) we obtain ct{Mw) = —0.195 , Cg{Mw) = 2.056 and Cio{Mw) = —4.415 (with 
Cq{M]v) defined according to the prescription given in |^). The coefficient Cg-^^ appearing 
in ( p.5|) includes the s dependence induced by the one-loop matrix element of the four-quark 
operators, and reads 



C9 



cg{mb) + [3ci(mfe) + C2{mb)] g{rhc, s). 



(3.7) 



where rhc = rric/mb and s = s/m^. In our numerical analysis, we use the expression for 
glthcs) given in [^]. An additional contribution to this decay mode comes from the long 
distance effects associated with on-shell and off-shell cc resonances. There are six known 
resonances that can contribute. They generate an additional term which has the same 
structure as O9. Hence, it is convenient to include the resonance contributions directly into 



by making the replacement 



3n 

g{mc, s) ^ g{mc, s) = g{mc, s) ^ 

rv — 

" v=.j/i),ii)',.. 



KvMvt{V 1+1-) 
s - M^, + iMvtv 



(3.8) 



where Ty and TiV l^l ) are respectively the total and partial decay width of the vector 
meson resonances normalized to ntb , while My is the normalized vector meson mass. The 
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values of the masses, widths and leptonic branching ratios of the six cc resonances can 
be found in [^], while the phenomenological factor ~ 2.3 is calculated by fitting the 
B J/ipK* amplitude to the experimental rate The replacement ( p.8| ) to model the 



long distance contributions from off-shell resonances is not a rigorous procedure |3^. Here 
we adopt this simple prescription in order to compare the long distance effects with the new 
physics short distance contributions. 

The low energy coefficients of the new operators in ( |3.2D can be determined in the same 
way as the SM coefficients. Due to the vectorlike nature of QCD, the matrix of anomalous 
dimensions for the set {O^} is the same as for the standard basis. Moreover, in the models 
that we will study in the next section, operators induced by squark exchange at the scale 
rriq do not renormalize. Therefore, to take into account the QCD effects on the new physics 
operators, the only new scale that needs to be introduced is m = rrij . We fix m = 100 GeV. 
We note that in this case the renormalization group evolution is still controlled by the 
QCD /5-function with five active flavors. Modifications to account for the case mj > rrit 
are straightforward. In the models discussed in the next section, C[{m) = 0. Since new 
physics generates Oj only at the loop level, we will also set Cj{rh) = and C^{mb) arises 
only from operator mixing. Then the set of equations that determines the low energy 
coefficients C[{mij) is the same than ( |3.(j| ) with the replacements c^{My/) — > Cjijh) = and 
1] ^ 7]' = as{rh) / asirrib) . In particular, we have included in C^Jijni,) also the additional long 
distance contributions induced by O'^ and in the same way as for cj^ in (|3.7|) . 

Finally, the evolution of the coefficients of the scalar operators in ( |3.3| ) is controlled by 
the anomalous dimensions of Of , Ofg , Of and Ofg : 70s = —4 . 



B. Inclusive rate and various observables. 



Neglecting non-perturbative (~ ^/fnl) corrections ||32| the inclusive h ^ Xgl'^ I decay 
width as a function of the invariant mass of the lepton pair = mf+i- is given by 



ds 3847r^ ^ ' ' V 



(3.9) 



where s = q'^/ml, rhi = mi/mb , A(a, b, c) = a? + h'^ + — 2{ab + bc + ac) and 



=4 1 + 



2m? 



i + |0^|2) Fi(S, m^) + ^Re {C;C',) F^is, ml 



+ 3 Re (07*09 + C'j*C'g) F^is, mj) - 3 Re (207*0^ + 20*0^ + 0*0^ 
+ (IO9P + |0^p + lOioP + \C[,\') F,{s,mlmf) - 12m,s (^1 - Re(0*oO 



V 1 2 



+ 6mt i\C,\' + - lOiol^ - \C\,\']F5is,m 



■I 1 2 



-"lOl 



3. 

+ 2^ 



1 - 



4m, 



O 



S|2 I \(jiS\2 



5 1 2 I \(jiS\2 



10 I 



10 I 



F.is.rr?^ 
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+ Q STTls 



1 - 



4m 



-"10 "-^lO 



s 

+ 6miRe(C(f Cio + C: 
- 6 Re(Cif C^q + Cfo* Cio) (^5(5 

The functions Fi read 



i)Re(crCf)+Re(cfo*C(i 
(F5(S,m2)_2m2 



(3.10) 



F4(s,m^,mf) 



1 - m^)2 - s(l + 14m^ + m^) 



s(4 + Ami + -5 



2(1 + m 
2(1 -m 

(1 — m^)^ — s{l + m 



,2\2 



;i - my + s(l + ) - 2r + A(l, s, mf) 



2mf 



(3.11) 



The forward-backward asymmetry A^-q is defined with respect to the angular variable 
ce = cos 6 where 6 is the angle of the lepton /~ with respect to the 6-direction in the 
center- of- mass system: 



Afb{s) 
We obtain : 



1 



dT{s)/ds 



dcf) 



ds dci 



■Tis.ce) 



d^ 

dc8———T{s,C0) 
^1 dsdcf) 



(3.12) 



Apb(5) = 3A1/2(1, 



where 



A(s) = sReiC;Cio-C'a*C- 



'10 



Amf A(s) 

s ns) 



+ 2(1 + ml)Re{c;Ci(, - C'^*C[^ 



(3.13) 



- 4m,Re(C7*C(o - Cl^C'^) + m Re [(cf + Cf) (2C; + 
+ nil Re [(cf * + m, Cf ) (2C7 + Cg)] . (3.14) 

Recently, polarization measurements of final state r leptons in the decay 6 — >■ r"*" r~ 
have been proposed as a useful tool in discerning physics beyond the SM p4| , p5| . Experimen- 
tal studies of the r polarization have been carried out by the four LEP collaborations, and 
from the analysis of the distributions of the r decay products in different decay modes, the 
r polarization was determined with an error of about 10% ||3^. For S-factory experiments 
the task of reconstructing the r polarization appears more challenging. For example, in 
contrast to the LEP environment, the energy of the decaying r is not a-priori known. How- 
ever, if a comparable sensitivity can be reached, then a measurement of the r polarization 
components could provide crucial informations for distinguishing among different sources of 
new physics. Following Ref. |35| we define the inclusive lepton polarization components by 
introducing three orthogonal unit vectors 



eL 



P 



IP- 



p. X P 

IPs X P- 



(3.15) 
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where p_ and are the three-momenta of /~ and of the s quark in the cm. frame of the 
/"'"/" system. The differential decay rate for b ^ Xg for any given spin direction n of 
the lepton /~ (n being a unit vector in the /~ rest frame) can be written as 



dr(s,n) _ 1 (dT{s) 
~ 2 



ds 



ds 



unpol 



1 + (Pl(s) eL + Pn(s) Gn + Pt(s) Gt) ■ n 



(3.16) 



where Pl, and Vt give the longitudinal, normal and transverse components of the 
/~ polarization as a function of s. The polarization asymmetries Vi{s) {i = L, N, T) are 
obtained by evaluating 

dr{ei, s)/ds — (ir(— ej, s)/ds 



Vi{s) 

dT{ei, s)/ds + dT{—ei, s) /ds 
for e, = bl , eN , ex . For the general form of the interaction ( p.4| ) we obtain 

1 I 4m? 



(3.17) 



'1- 



12Re(C;Cio + C;*Cg (1- 



m 



2\2 



S(sj V s 

+ 2 Re (C*Cio + C';C[q) [(1 - mlY + s{l + m^) - 2s^ 
- 12 rriss Re (2 C;C[o + 2 C!;Cw + C;C[o + C'^Cw) 



s{l + 



- 3s 



(1 - s + ml) Re {cfC^, + CfC[f) - 2m, Re {cfC[l + CfC- 



6 rhi m, (1 + s - mi) Re (^9^*^10 + C 



'10 



(3.18) 



Pn(s) 



2lm{-C;C[o + C'*Cw) rhs 



+ Im(C;Cio - C^7*C^io) (1 + + 2 m; yilm (C*Cio - C';C[ 



'lOJ 



VI Im (C'f (2 C7 + Cg) + Cf (2 + C^) + (C;f Cio + C7fo*C(o 



m.. 



Im (Cf (2 C7 + C9) + * (2 + C^; 



•-^10 "-^lO 



^5 *^ 



, (3.19) 



Pt(S) 



1 37r 



-|(|C^7r-|C^;P) (l-m^)^ 



S(s)2v^ 

+ Re (C*oC9 + CfoC^) (1 - m^) - 8 m,Re (C;*C9 - C^C'^ -s[\C^\^- \C'^\' 
- 4 Re {C;C^ - C'*C'^) (1 + m^) + 2 Re (C^Cio + C;*C(o) (1 - ml) 



s 

+ 2 



4mf ' 



Re (Cf *Cio + Cf C(o) + Re (cf Cio + *C( 
+ Re (C^f (2 Cr + Cg) + Cfo* (2 C, + C^)) 
+ m. Re (Cfo* (2 C7 + Cg) + Cjf (2 + C'^ 



(3.20) 
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The possibility of measuring the longitudinal polarization of the r in the decay 
b XgT^ T~ was first proposed in [34|. This analysis was extended to include the other 
two polarization components Pn and Vt in p5[. The SM results are recovered from our 



formulae by setting the C- and the Cg [q coefficients to zero. In this limit our expressions for 
and Pn(s) agree with the results given by Kriiger and Sehgal and when we set 
rris = in (|3.18| ) we find agreement with Pl as given in [^]. However, the SM limit of Vt 
disagrees with eq. (5.5) in ||3^ for the factor of two multiplying the term Iie{CjCiQ) in the 
third line of QOOl) . 

To derive numerical results for the inclusive branching ratios and for the averaged values 
of the asymmetries, we minimize long distance effects by imposing cuts on the dilepton 
invariant mass. For the muon channel we select the region below the resonances s < 0.4 , 
while for the tau channel we require s > 0.6 , above the ip' . The remaining effects of the 
four additional cc resonances in the tail of the invariant mass distribution for b XsT~^ t~ 
are not very large. 

At leading order, the SM results for the various quantities are as follows. The SM 
inclusive branching ratio for muons is predicted to be BR(6 — > yU+yU^)^-<o.4 = 4.3 x 10^^ . 
The purely short distance contribution yields in the same region H )s<0.4 — 

3.9 X 10~^ . Due to phase space suppression and to the different cut, the tau production rate 
is more than one order of magnitude smaller. We find BR(6 — > X^ r"'"r~)s->o.6 = 1-5 x 10"'' 
and a similar number for the purely short distance contribution. Below our cut, the average 
value of the yU forward-backward asymmetry is rather small (^fb)j<0 4 ~ —0.01. This is 
due to the fact that the asymmetry changes sign in this region, and in taking the average 
large cancellations occur. For the r lepton the asymmetry is larger (^fb)5>0 6 ~ —0.13. 
Only the longitudinal polarization asymmetry Pl is significant for the : (Pl)s<0 4 ~ —0.57 
while all the three components are sizeable for r . We find the following average values 
(^£).>o.a = -0-34, = -0.40, = 0.05. 

Notice that the normal component Pn, though small, is a T-odd quantity and it is 
considerably larger than the corresponding normal polarization of leptons in vr^/i"// 



or —>■ TT^fi^fi . However, this should not be interpreted as signaling CP violation. In 



fact, since we assume real couplings, the leading contribution comes from the absorptive part 
of the effective coupling Cg (and Cg in the new physics analysis), which are dominated by 
the cc real intermediate states (cf. ( ^.8D ). It is in principle possible to remove this background 
to CP violating effects by measuring the difference between the Pn components for example 
in and B~ decays. This is because CP violating phases yield asymmetries with the same 



sign, while phases originating from strong interactions cancel off |38 



IV. SUSY WITHOUT R-PARITY 



In this section we present a short introduction to SUSY models where R-parity is not 
imposed, and L is (mildly) violated already at the renormalizable level. Next we will embed 
these models in the framework of an Abelian horizontal symmetry that will provide 
us with the additional theoretical constraints needed to derive numerical predictions. As 
it was discussed in 0, this procedure has the additional advantage of ensuring that in the 
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particular models we will discuss, the present constraints on the R-parity violating couplings 
are satisfied. 



A. R-parity violating couplings 

The field content of the SM together with the requirement of SU{2)l x U{1)y gauge 
invariance, implies that at the renormalizable level the most general Lagrangian possesses 
additional accidental U{1) symmetries, corresponding to conserved baryon and lepton fiavor 
{Li) quantum numbers. The conservation of i? , Lj and hence of total lepton number 
(L = naturally explains nucleon stability as well as the non observation of L and Lj 

violating transitions. In SUSY extensions of the SM, additional gauge and Lorentz invariant 
terms are allowed, which violate B, Li and L . Denoting collectively by Ha {ot = 0,1,2,3) 
the supermultiplets containing the down-type Higgs and the left-handed lepton doublets, 
which transform in the same way under the gauge group, the following Li and L violating 
superpotential terms arise 

Wf; = l^aHaHu + Xa/3k Ha Hp 1% + \' ajk Ha Qj . (4.1) 

Here Qi and df denote the quark doublet and down-quark singlet superfields, are the lepton 
singlets and H^ contains the up-type Higgs field. There are also renormalizable terms which 
violate B, = X'-ji^Ui dj dl , and physics at some large scale M\ can induce additional 

dimension 5 B and L violating terms like (F^j^-^j/MA) Ha QiQjQk + ■ ■ ■ ■ 

To forbid the dangerous dimension 4 terms, a parity quantum number R = (_i^3B+l+2S 
{S being the spin) is assigned to each component field, and invariance under R transforma- 
tions is imposed. However, even if suppressed by the Plank mass, the R-parity conserving 
dimension 5 terms can still induce too fast proton decay unless V < 10~^ etc. From a phe- 
nomenological point of view, the first priority is to ensure the absence of operators leading 
to fast nucleon decay, and in this respect other discrete symmetries can be more effective 
than R ^0|. These interesting alternatives forbid dimension 4 and 5 B violating terms but 
do not imply the same for the L non- conserving terms. Since a mild violation of L can 
be phenomenologically tolerated, SUSY extensions of the SM with highly suppressed B vi- 
olation but without R parity and without L number, represent interesting alternatives to 
the Minimal Supersymmetric Standard Model (MSSM). We henceforth assume that B is 
effectively conserved, and we concentrate only on the L- violating terms contained in (|4.1|). 



The first term in (4J.) can mix the fermions with the Higgsinos, resulting in too large 
neutrino masses P,^T|-^]. In Ref. |^ this problem was solved by assuming that the down- 
type Higgs transforms differently from the lepton doublets under a horizontal symmetry. The 
different charge assignments can generate enough suppression of the lepton mixing with Hu ■ 
However, this does not occur in the models discussed in the next section (and in Ref. 0]). 
In fact in these models the down-type Higgs and the r doublets have the same horizontal 
charge, and in the model labeled below as Model II, this is true also for the muon doublet. 
We can still avoid generating too large neutrino masses under the assumption that the soft 
SUSY breaking terms are universal. As it is discussed in [^], this assumption implies the 
following: 

(i) Modulo small violation of the universality conditions at the electroweak scale induced 
by the renormalization group running of the soft SUSY breaking parameters, the 
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combination Hd = fiaHa/ a/aw^ corresponds to the down-type Higgs. Namely its 
scalar component is the only other field, besides H^, acquiring a non- vanishing vacuum 
expectation value. 

{ii) In the basis {Ha, Li} (where Li with z = 1, 2, 3 denote the three combinations orthog- 
onal to Hd) the couplings Xoij and A'oy are respectively the Yukawa couplings Y^j of 
the leptons and Y^^ of the down-type quarks. 
As it will become clear in the next subsection, this short discussion is relevant only for the 
subset of the multiplets that carry the same minimum charge. For doublets having hor- 
izontal charge assignments different from the minimum one (as for example for the electron 
doublet) the horizontal symmetry by itself defines to a very good approximation the physical 
fields. 

Once the fields are rotated to the physical basis, the L violating trilinear terms contained 
in (|4.1|) read 

Kjk Li Lj + X'ijk Li Qj , (4.2) 

where Xijk = —Xjik due to the antisymmetry in the SU{2) indices. Several of the A and A' 
couplings are strongly constrained by the existing phenomenology The best limits are 
for couplings involving fermions of the first two generations (i, j, fc = 1, 2) while for couplings 
involving more than a single third generation field the existing limits are much weaker, and 
in some cases no bounds exist to date. This situation is interesting since in general models 
that can explain the observed fermion mass hierarchy also predict that R-parity violating 
couplings involving third generation fields are the largest ones. 



B. R-parity violation in the framework of horizontal symmetries 



In order to evaluate the effects of the new R-parity violating interactions, we need to 
estimate quantitatively the coefficients A and A'. We work in the framework of the super- 
symmetric models with horizontal symmetries that have been thoroughly investigated in 
|39| . These models successfully predict the order of magnitude of the fermion masses and 
CKM angles, and can also explain the suppression of L violation 0] and B violation in 
SUSY models without R-parity. Assuming a horizontal U{1) symmetry allows us to esti- 
mate the size of the L violating couplings, and to work out numerical predictions for various 
observables measurable in B decays. 

In the models we are interested in, there are no additional fields in the low energy 
spectrum with respect to minimal SUSY. However, a charge H{%1)) of an Abelian horizontal 
symmetry Ti = U {1) h '^^ assigned to each supermultiplet ip . Tiis explicitly broken by a small 
parameter e with charge H{e) = —1 giving rise to a set of selection rules for the effective 



couplings of the low energy Lagrangian [^. If we assume that each of the lepton, quark and 
Higgs superfields carries positive or zero charge, the selection rule relevant for the present 
discussion is that the effective coupling gabc for a general trilinear superpotential term ipai^bi^c 
is of order gabc ~ £^i'^'^)+^i'^b)+H{ipc) _ Therefore the leptons and down-type quarks Yukawa 
couphngs are respectively of order ~ ^Hi^,)+H(L^)+H(i^^} yf. ~ ^m^)+Hm+H{d'^) 
(rotation to the exact quark mass eigenstate basis does not affect these order of magnitude 



14 



estimates 0). Most of the L- violating couplings in ( [4.2| ) are further suppressed with respect 
to the corresponding Yukawa couplings. They can be estimated as 



1 /2 

\ COS l) j 



and 



A' ~ yd H{U)-H{^^) _ / 2v^2G£\ ' H{d^^)-Hid^^)+H{L,)-Hi'f^) ^-^ 

These equations show that 
{i) the couplings A and A' involving fermions of the third generation 

are respectively enhanced by and rrib ; 
(ii) like the lepton and down quark Yukawa couplings, the A, A' couplings increase with 
tan P . 

In order to give a numerical estimate of the couplings, we need a set of charges and a value 



for the ?i-symmetry breaking parameter e . In the model discussed in |]39|, e ~ 0.22 is fixed 
by the magnitude of the Cabibbo angle, while the quark, lepton and Higgs charges are chosen 
to reproduce the values of the fermion masses and CKM mixing angles. Besides reproducing 



the measured values, the model has some predictivity in the quark sector it yields 



estimates for ratios of neutrino masses and most important in the present context, it 
ensures that the L- violating couplings in ( [4.2| ), (^73|) and ( |4.4| ) are safely suppressed below 



the present experimental limits The following 7i-charge assignments [^] fit the order of 
magnitude of all the quark masses and CKM mixing angles 

Qi Q2 Q3 dl dl dl ul ul ul 
(3)(2)(0) (3) (2) (2) (3)(1)(0) (0) (0) . (4.5) 

For the leptons, we will use two different sets of charges that fit well the order of magnitude 
of the charged lepton masses. We also use a different value for the squark masses for 
each set : 

Li L2 L3 il il m-iiGeY) m^- (GeV) 

Model I : (4) (2) (0) (4) (3) (3) 100 170 

Model II : (3) (0) (0) (5) (5) (3) 100 350 . (4.6) 



The charge assignments and the sfermion mass values listed in (|4.5| ) and ( [4.6|) , together with 
e ~ 0.22 completely define the two models and allow us to estimate the order of magnitude 
of the various decay rates. The charges of Model I coincide with the charges of the "master 
model" of [Q. While in this model new physics effects are induced dominantly by the new 
operators {O'} in ( |3.2| ) arising from squark exchange, in Model II the leading effects are due 
to the scalar operators {O^} in (|3.3| ) induced by slepton exchange. The choices (|4.5| ) and 
( [4.6[ ) for the horizontal charges are not unique. Since the Yukawa interactions are invariant 
under a set of U{1) symmetries such as hypercharge or lepton number, it is always possible 
to shift the if-charges of any amount proportional to one of the corresponding f/(l) quantum 
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numbers, without affecting the predictions for the masses and mixing angles.^ In particular, 
a shift proportional to L : H{Li) H{Li) + n, H{if) H{if) - n and Hi^p) H{^p) for 
all the other fields, has the effect of suppressing (for n > ) all the L violating couplings 
in ( [4.2|) , ( [4 .31) and ( [4 .41) by a factor of . We have found that already for n = 1 the 
couplings are enough suppressed so that no signal of new physics can be detected in the 
experimental quantities we are considering in this paper. Notice also that Model II can be 
derived from Model I by means of shifts proportional to lepton flavor numbers: rig = — 1 , 
= —2 , rir = . This has the effect of enhancing some of the A couplings without affecting 
the charged lepton masses. Of course, in Model II the predictions for neutrino masses and 
mixings will differ from the predictions of Model I M . 



C. Coefficients and rates for the various decays 



In SUSY models without R-parity, h r decays can proceed through the exchange of 
sleptons and/or squarks, yielding significant enhancements over the SM rates. In this section 
we give the expressions for the various coefficients of the effective operators contributing to 
the decays ^ t u, B ^ , b XgU u and 6 — > r+ including the new physics 

contributions. 

For the decay B~ — > r z/ the coefficients appearing in (p.3|) and (p.4|) read 



^ 31fc^i3fc 

AGf m?, 
2Gf mf 



V2V, 



ub 



(4.7) 



where a sum over the repeated index k is left understood. The index i refers to the final 
state neutrino flavor which, as already said, can be different from z/^. 
For the decay Bq — > r+ the coefficients in ( p.6| ) and ( |2.?| ) read 



CI 



CI' 



C\ 



A:33-^fc3(j + 



AGprn^- 



AGprnf 

Ik 



^ 3fc3^3fc(? 

'sGf 



+ 



CI 



Qk 



SM 



(4.8) 



where 



CI 



SM 



IS given m 



For the decay b —>■ XqUiUj the coefficients appearing in ( 2.10 ) and ( |2.11 ) read 



•'■In models where the 7i symmetry is gauged, shifts of this kind cannot be arbitrary but must 
respect the constraints from anomaly cancellation. 
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2GFmj 3 

^ ik3^jkq 



CI 



SM 



dk 



(4.9) 



where 



SM 



IS given m (pl2|) . Finally, for the decay 6 X, r+ r , the coefficients of the 



new operators {O'} {i = 2, 7, 9, 10) in ( p.2| ) and {O } in ( |3.3| ) at the scale m 

^'k23^'k22 



^ 3fc3^ 3fc2 

4GFm?^ 



m^- are 







[m) = 






[m) = 


G',{m) = 




[m) = 


Cl{m) = 


•-^10 


[m) = 


Cf (m) = 


r^s ' 


'jfi) = 



2Gf 



fk 



^'*k23^k33 

2GFml 



(M3) 
(M3) 



(4.10) 



As regards the coefficients of the standard operators in ( |3.1| ), they are affected by the new 
physics only at the loop level. Since these are subleading effects we will neglect them. 

Few comments are in order. The antisymmetry in the ffist two indices of the A couplings 
forces /c 7^ 3 in the coefficients of the scalar operators. In Model I, because of the lepton 
charge assignments (|4.6| ) this results in a strong suppression of the scalar couplings. We 
have for example Gg ~ Gg so that the leading new physics effects are due to the operators 
induced by squark exchange. In Model II, being H{L2) = H{L^) = there is no suppression 
of the scalar couplings from the ?i-symmetry. Then the large value of the squark masses 
rriq = 350 GeV suppresses Gg and C(q down to ~ eGg so that in this model slepton exchange 
gives the dominant effects. As regards the b ^ Xg /i~ decay channel, the charge difference 
H[L2) — H{L'i) = 2 of Model I implies a strong suppression of the A'A' couplings. In contrast, 
for the couplings involving the combination A'A the antisymmetry now allows = 3 , with 
the noticeable result that the scalar couplings in the /i channel are enhanced with respect to 
the r channel. Still the enhancement is only of about a factor e"^ so that we cannot expect 
particularly large new physics effects. In Model II there is no suppression of the squark 
exchange operators for b Xg fi'^ fi~ with respect to the r channel. Thus we can expect 
that Gg and G[q will still produce observable signals of new physics. Finally, C2 that controls 
the relevance of the new long distance contributions, is rather small in both models. This is 
due to the charge H{Q2) = 2 which yields a relative factor 5^ {mq/miY ~ 0.001 (Model I)- 
0.03 (Model II) with respect to the Gg short distance contribution. Therefore we can expect 
that the resonance peaks will be smeared off by the dominant new physics short distance 
effects. 
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V. RESULTS 



The main results of our analysis consist of a set of numerical predictions for several 
observables in the presence of new physics from SUSY models without R-parity. We recall 
that the theoretical framework we adopted is a straightforward extension of successful models 
for fermion masses and CKM mixing angles, and gives predictions for the R-parity breaking 
couplings which are consistent with all the present experimental constraints 0. 

We have studied several observables measurable in 6 — r , z/^ and 6 — /x transitions 
within two different models, that were defined in ([4.5| ) and (|4.6| ). In Model I, the horizontal 
charges coincide with the charges of the "master model" discussed in The sfermions 

masses are rrii = 100 GeV and rriq = 170 GeV. In this model the dominant new physics effects 
come from squark exchange which generate the leading new effective operators. Model II was 
introduced in order to study the effects of sleptons. It is defined in terms of a different set of 
lepton horizontal charges {WMj- The sfermion masses are mj = 100 GeV and rriq = 350 GeV. 

Varying the values of the new physics parameters results in the following scaling behaviors 
of the branching ratios: 

(i) In both models, the branching ratios scale as (tan/3)^ and as e^", where n is an 
arbitrary shift of the if-charges of the leptons proportional to lepton number (see the 
discussion at the end of section IV-B). Our results are given for n = and tan/3 = 1 . 
{ii) In Model I the branching ratio for b XsT'^ t~ scales as (170 GeV/rriq)^. Hence, for 
light squarks {rriq ~ 100 GeV) and moderate values of tan/5 ( > 2) Model I can predict 
values of the branching ratio up to few xlO~^ . In Model II the rate for b ^ XgT^ t~ 
scales as (100 GeV/mj)^ so that the main enhancement with respect to our predictions 
(see table I) can only come from tan (3 > 1. 
Our results are collected in tables I and II, and in figures 1-8. 

Table I lists the predictions for decays involving the transition 6 — r , i/^. The first five 
lines give the results for the decays ^ t u, Bq ^ r+ and b X u u, while the results 
for the branching ratio, forward-backward asymmetry and longitudinal polarization asym- 
metries in b ^ Xg T~ are given in the remaining entries. Table II collects the results for 
the corresponding processes with final state muons. In both tables, the first column lists the 
SM predictions for the various observables (computed in the leading order approximation) 
the second column lists the predictions of Model I while the results for Model II are listed 
in the third column. 

Since the two decays b ^ Xg fj.~^ fi~ and 6 — > r+ are affected by large long distance 
effects, to single out the short distance contributions we have applied cuts on the dilep- 
ton invariant mass. We study b ^ X^t'^ t~ in the region above the ip' (■§ > 0.6) while 
b Xs/i+/i~ is analyzed below the resonance region ( s < 0.4). A comparison between 
the total inclusive branching ratios BR(6 — Xg /~)nocut , the branching ratio in the region 
within the cuts and the kinematic limits BR(6 —>■ Xgl'^ l~)s<oA(s>o.a) and the short distance 
contribution in the same region BR'''^(fe — > Xs l~) shows the effects of the cuts on the total 
rates, and their effectiveness in isolating the interesting contributions. 

From the results in table I, it is apparent that in most cases the decays B^ ^ tv and 
b —>■ Xq V V are not very sensitive to the sources of new physics we are analyzing here. In both 
our models these decays do not show any significant enhancement with respect to the SM 
rates. New physics from Model I can enhance the rates for Bg — >■ t~ and B^ t^t~. 
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However, the rates remain very small, and because of the large theoretical uncertainty 
related to fs it is not obvious that these signals could be unambiguously identified. In 
Model II the branching ratio for Bg — > r"*" r~ increases by two orders of magnitudes, up to 
~ 10~^. This process cannot be searched for at 5-factories running at the T(45') . However, 
hadron colliders might be able to detect this enhancement, depending on the efficiency in 
identifying the r's. In this model we also observe a two orders of magnitude enhancement for 
Bd T~ which appears more promising for new physics searches at future B-factories. 

In the muon channel, the corresponding decays Bg — > /x+ /x" are also sensitive to new 
physics effects from both models. However, even if the decay rates are enhanced by two 
orders of magnitude, the branching ratios are still only at the level of ~ 7 x 10~^ for Bg and 
~ 3 X 10-^ for Bd. 

As regards the decay b ^ XgfJ^^ we see that no signal of new physics is expected 
in Model I. The rate, the forward-backward asymmetry and the longitudinal polarization 
asymmetry remain at their SM values. This decay is somewhat more sensitive to new 
physics from Model II, which induces a factor of two enhancement of the short distance 
contributions. However, as it clear also by inspecting figs. 6-8, below the cut s < 0.4 the 
overall effects are very likely too small to be unambiguously identified above the theoretical 
and experimental uncertainties. 

We turn now to the decay b —>■ XgT~^ t~ . In both our models the branching ratio for this 
decay is enhanced by more than one order of magnitude, at a level that could be observable 
with a good r identification efficiency. Notice that with our choice of new physics parameters. 
Model I and Model II both predict very similar rates (see also fig. 1) even if the respective 
enhancements are induced by effective operators of quite a different nature. Therefore, 
even if such a large signal of new physics will be observed, it would not be possible to 
identify which kind of new physics is producing the effect just from a measurement of the 
decay rate. In contrast, we see that the forward-backward asymmetry A^b (fig- 2) and the 
longitudinal polarization asymmetry Vh (fig- 3) could provide the additional information 
needed to disentangle the different effects. 

In Model I Ap-g is almost at the 20% level, and opposite in sign with respect to the SM 
in the whole kinematic region. In contrast, in Model II this asymmetry is vanishingly small. 
On the other hand in Model II the longitudinal component Pl is rather large, about twice 
the SM prediction, while in Model I it remains close to the SM value. 

As regards the other two polarization components, in both models the transverse asym- 
metry Vt is about a factor of three smaller than in the SM. The T-odd component which 
in the SM is at the 5% level, is practically zero in both our models. This can be traced back 
to the fact that the new (real) short distance contributions dominate over the absorptive 
part of the decay rate related to on shell cc intermediate states. 

VI. CONCLUSIONS 

In this paper we have studied rare b decays into leptons of the third generation. These 
decay modes are well suited to study sources of new physics which couple more strongly 
to third family. We have first discussed a general framework for studying effects beyond 
the SM, and we have introduced general four-fermion amplitudes for the decays B~ — > r P, 
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Bq — > r"*" T~ and b ^ XqU u. We have also defined the effective Hamiltonian for the decay 
b Xg r+ in terms of an enlarged operator basis. 

We have applied our results to the study of SUSY models without R-parity and without 
L number. In these models new contributions to the decays appear already at the tree 
level, through new effective operators generated by squark and slepton exchange having 
a different structure than the SM ones. In order to derive numerical predictions for the 
various observablcs, we have embedded SUSY without R-parity in the framework of models 
for fermion masses based on Abelian horizontal symmetries. This allowed us to estimate the 
order of magnitude of the various R-parity violating couplings. 

We have carried out a numerical study of two representative models, in which new physics 
effects arise from two different sets of effective operators, induced respectively by squarks 
and by sleptons exchange. We found that the most sensitive among the processes involving 
the b —>■ T transition are the decays Bg —>■ t~ {q = d,s) and b —>■ XgT^ t~ , that can 
be enhanced up to two orders of magnitude over the SM rates. If such an enhancement 
is observed, additional measurements of the forward-backward asymmetry and of the r 
longitudinal polarization in b ^ Xg r~ can be very helpful in identifying the kind of 
underlying new physics. We have confronted the predictions for decays into final state 
taus with the corresponding results for decays into muons. The decay 6 — > z^"*" /x~ is not 
sensitive to this kind of new physics, and we found that only the Bg /x+ /x~ decay modes 
show enhancements comparable to the r channel. However, even with the new physics 
contributions, the overall rates for these decays remain rather small. 
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APPENDIX A: INPUT PARAMETERS 

rub = 4.8 GeV, rUc = 1.4 GeV, = 0.2 GeV, rrit = 176 GeV , 
= 0.106 GeV, = 1.777 GeV, Mw = 80.2 GeV, 

Vt, = L Vts = Vcb = -0.040, = 0.009, Kb = 0.003, 
fe^ = 200 MeV, /s, = 230 MeV, = 5.3 GeV, tb = 1.6 ps 
Aqcd = 225 McV, fi = mh, a{Mz) = 1/129, sin^ 9^ = 0.23, 
BR{B^XJui) = 10.4%. 
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TABLE I. Predictions for the various decay rates and asymmetries for b ^ t transitions in the 
standard model and in the R-parity violating models discussed in the text. Model I is sensitive 
to operators generated by squark exchange. The lepton horizontal charges are H{L) = (4,2,0), 
Hii") = (4,3,3) while the SUSY masses are m^- = 100 GeV and rug = 170 GeV. Model II is 
sensitive to operators generated by slepton exchange, with horizontal charges H{L) = (3,0,0), 
i7(r) = (5,5,3) and SUSY masses = 100 GeV and ruq = 350 GeV. In both models the value 
of the horizontal symmetry breaking parameter is e = 0.22 . 



Process 


Standard Model 


Model 1 




Model 2 


BR {B- ^ T- v) 
BR {Bs r+r") 
BR (Bd r+r") 
BR{b^ Xsui?) 
BR {b ^ Xdiyu) 


7.1 X 10-^ 
9.1 X 10"^ 

4.3 X 10"^ 

4.4 X 10"^ 
2.7 X 10-6 


7.2 X 10-5 
5.7 X 10-6 
1.9 X lO-'^ 
6.7 X 10-5 
3.9 X 10-6 




7.4 X 10-5 
1.8 X 10-^ 
6.3 X 10-6 
5.0 X 10-5 
3.0 X 10-6 


BR {b^XsT+T-U,^, 

BR (6^X,r+T-),>o.6 
BR^<l(^^^^^+^-),^^^ 

(^Fb)s>o.6 

\ ' L / i>0.6 
\ ' T /s>0.6 
( "^N )|>0.6 


4.9 X 10-6 

1.5 X 10-^ 

1.6 X lO"'^ 
-0.13 
-0.34 
-0.40 

0.05 


9.6 X 10-6 
4.1 X 10-6 
4.1 X 10-6 
0.18 
-0.40 
-0.13 
0.00 




1.0 X 10-5 
4.6 X 10-6 
4.6 X 10-6 
-0.03 
-0.68 
-0.14 
0.01 


TABLE II. Predictions for the various decay rates and asymmetries for b ^ fi transitions in 
the standard model and in the R-parity violating models discussed in the text. Model I is sensitive 
to operators generated by squark exchange. The lepton horizontal charges are H{L) = (4,2,0), 
H{i'') = (4,3,3) while the SUSY masses are rrif = 100 GeV and rriq = 170 GeV. Model II is 
sensitive to operators generated by slepton exchange, with horizontal charges H{L) = (3,0,0), 
H{1'^) = (5,5,3) and SUSY masses m-^ = 100 GeV and rriq = 350 GeV. In both models the value 
of the horizontal symmetry breaking parameter is e = 0.22 . 


Process 


Standard Model 


Model 1 




Mcxlel 2 


BR {B- ^ ij.- u) 

BR (B, ^ fj+ii-) 
BR {Bd ^ /X+/X-) 


3.2 X lO"'^ 

4.3 X 10"^ 
2.1 X 10~^° 


3.2 X 10-^ 
7.9 X 10-^ 
2.9 X 10-^ 




3.3 X 10-'^ 
7.2 X lO-'^ 
2.7 X 10-^ 


BR (6 Xs^+^")n„^„t 
BR (6 Xs n^n~)s<oA 
BR^d(6^X,^V")i<o.4 

(^Fb)j<o.4 

\ ' L / s<0.4 


3.1 X 10"*^ 
4.3 X 10^6 
3.9 X 10-6 
-0.01 
-0.57 


3.1 X 10-^ 
4.5 X 10-6 
4.1 X 10-6 
0.00 
-0.56 




3.4 X lO-'' 
8.3 X 10-6 
7.7 X 10-6 
0.08 
-0.73 
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FIG. 1. Predictions for the differential branching ratio BR(6 XgT'^ t~) as a function of s 
in the standard model (solid), in Model I (dashed) and in Model II (dash-dotted) discussed in 
the text, including the long distance contribution. Model I is sensitive to operators generated by 
squark exchange. The lepton horizontal charges are H{L) = (4,2,0), H{1'^) = (4,3,3) while the 
SUSY masses are mj = 100 GeV and ruq = 170 GeV. Model II is sensitive to operators generated 
by slepton exchange, with horizontal charges H{L) = (3,0,0), H{F) = (5,5,3) and SUSY masses 
mi = 100 GeV and rriq = 350 GeV. In both models the value of the horizontal symmetry breaking 
parameter is e = 0.22 . 
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FIG. 2. Predictions for the forward-backward asymmetry ApB for the r lepton as a function of 
s in the standard model (solid), in Model I (dashed) and in Model II (dash-dotted) discussed in 
the text. The new physics model parameters are as in fig. 1. 
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FIG. 3. Predictions for the longitudinal polarization Vh for the r lepton as a function of s in 
the standard model (solid), in Model I (dashed) and in Model II (dash-dotted) discussed in the 
text. The new physics model parameters are as in fig. 1. 
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FIG. 4. Predictions for the transverse polarization Vt for the r lepton as a function of s in the 
standard model (solid), in Model I (dashed) and in Model II (dash-dotted) discussed in the text. 
The new physics model parameters are as in fig. 1. 



0.10 



0.08 



0.06 



. 04 



0.02 








0.6 0.7 ^ 0.8 0.9 



FIG. 5. Predictions for the normal polarization Vn for the r lepton as a function of s in the 
standard model (solid), in Model I (dashed) and in Model II (dash-dotted) discussed in the text. 
The new physics model parameters are as in fig. 1. 
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FIG. 6. Predictions for the differential branching ratio BR(6 XgfJL^ as a function of s 
in the standard model (solid), in Model I (dashed) and in Model II (dash-dotted) discussed in the 
text, including the long distance contribution. The new physics model parameters are as in fig. 1. 
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FIG. 7. Predictions for the forward backward asymmetry Ayb for the /x lepton as a function 
of s in the standard model (solid) , in Model I (dashed) and in Model II (dash-dotted) discussed in 
the text. The new physics model parameters are as in fig. 1. 
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FIG. 8. Predictions for the longitudinal polarization for the fi lepton as function of s in the 
standard model (solid), in Model I (dashed) and in Model II (dash-dotted) discussed in the text. 
The new physics model parameters are as in fig. 1. 



30 



